Uncovering dark matter in cancer by identifying epigenetic drivers.

The complex relationship between chromatin accessibility, transcriptional regulation, and cancer transitions presents a daunting puzzle. Terekhanova et al. created a pan-cancer epigenetic and transcriptomic atlas at single-cell resolution, yielding important insights into the underlying chromatin architecture of cancer transitions and novel discoveries with the potential to advance precision medicine.

3D genomic features across >50 diverse cell types reveal insights into the genomic architecture of childhood obesity.

The prevalence of childhood obesity is increasing worldwide, along with the associated common comorbidities of type 2 diabetes and cardiovascular disease in later life. Motivated by evidence for a strong genetic component, our prior genome-wide association study (GWAS) efforts for childhood obesity revealed 19 independent signals for the trait; however, the mechanism of action of these loci remains to be elucidated. To molecularly characterize these childhood obesity loci we sought to determine the underlying causal variants and the corresponding effector genes within diverse cellular contexts. Integrating childhood obesity GWAS summary statistics with our existing 3D genomic datasets for 57 human cell types, consisting of high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq, and RNA-seq, we applied stratified LD score regression and calculated the proportion of genome-wide SNP heritability attributable to cell type-specific features, revealing pancreatic alpha cell enrichment as the most statistically significant. Subsequent chromatin contact-based fine-mapping was carried out for genome-wide significant childhood obesity loci and their linkage disequilibrium proxies to implicate effector genes, yielded the most abundant number of candidate variants and target genes at the BDNF, ADCY3, TMEM18 and FTO loci in skeletal muscle myotubes and the pancreatic beta-cell line, EndoC-BH1. One novel implicated effector gene, ALKAL2 - an inflammation-responsive gene in nerve nociceptors - was observed at the key TMEM18 locus across multiple immune cell types. Interestingly, this observation was also supported through colocalization analysis using expression quantitative trait loci (eQTL) derived from the Genotype-Tissue Expression (GTEx) dataset, supporting an inflammatory and neurologic component to the pathogenesis of childhood obesity. Our comprehensive appraisal of 3D genomic datasets generated in a myriad of different cell types provides genomic insights into pediatric obesity pathogenesis.

Master regulator activity QTL protocol to implicate regulatory pathways potentially mediating GWAS signals using eQTL data.

Here, we present a protocol to identify transcriptional regulators potentially mediating downstream biological effects of germline variants associated with complex traits of interest, which enables functional hypothesis generation independent of colocalizing expression quantitative trait loci (eQTLs). We describe steps for tissue-/cell-type-specific co-expression network modeling, expression regulator activity inference, and identification of representative phenotypic master regulators. Finally, we detail activity QTL and eQTL analyses. This protocol requires genotype, expression, and relevant covariables and phenotype data from existing eQTL datasets. For complete details on the use and execution of this protocol, please refer to Hoskins et al.1.

Relationship between ABO Blood Group Alleles and Pancreatic Cancer Is Modulated by Secretor (FUT2) Genotype, but Not Lewis Antigen (FUT3) Genotype.

BACKGROUND: In Western populations, pancreatic ductal adenocarcinoma (PDAC) risk has been found to be greater among individuals with non-O blood types than those with O blood type. However, the association has not been fully evaluated with respect to FUT2 (determining secretor status) and FUT3 (determining Lewis antigens) status, two biologically important genes in the expression of ABO blood groups with PDAC. METHODS: We examined interactions in data from 8,027 cases and 11,362 controls in large pancreatic cancer consortia (PanScan I-III and PanC4) by using genetic variants to predict ABO blood groups (rs505922 and rs8176746), secretor status (rs601338), and Lewis antigens (rs812936, rs28362459, and rs3894326). Multivariable logistic regression was used to estimate ORs and 95% confidence intervals (CI) of the risk of PDAC adjusted for age and sex. We examined multiplicative interactions of ABO with secretor status and Lewis antigens by considering each product term between ABO and secretor and between ABO and Lewis antigens individually. RESULTS: We found that the increased risk associated with non-O blood groups was somewhat stronger among secretors than nonsecretors [ORs, 1.28 (95% CI, 1.15-1.42) and 1.17 (95% CI, 1.03-1.32) respectively; Pinteraction = 0.002]. We did not find any interactions between ABO and Lewis antigens. CONCLUSIONS: Our large consortia data provide evidence of effect modification in the association between non-O blood type and pancreatic cancer risk by secretor status. IMPACT: Our results indicate that the association between ABO blood type and PDAC risk may vary by secretor status, but not by Lewis antigens.

Integrated Physiology of the Exocrine and Endocrine Compartments in Pancreatic Diseases: Workshop Proceedings.

The Integrated Physiology of the Exocrine and Endocrine Compartments in Pancreatic Diseases workshop was a 1.5-day scientific conference at the National Institutes of Health (Bethesda, MD) that engaged clinical and basic science investigators interested in diseases of the pancreas. This report provides a summary of the proceedings from the workshop. The goals of the workshop were to forge connections and identify gaps in knowledge that could guide future research directions. Presentations were segregated into six major theme areas, including 1) pancreas anatomy and physiology, 2) diabetes in the setting of exocrine disease, 3) metabolic influences on the exocrine pancreas, 4) genetic drivers of pancreatic diseases, 5) tools for integrated pancreatic analysis, and 6) implications of exocrine-endocrine cross talk. For each theme, multiple presentations were followed by panel discussions on specific topics relevant to each area of research; these are summarized here. Significantly, the discussions resulted in the identification of research gaps and opportunities for the field to address. In general, it was concluded that as a pancreas research community, we must more thoughtfully integrate our current knowledge of normal physiology as well as the disease mechanisms that underlie endocrine and exocrine disorders so that there is a better understanding of the interplay between these compartments.

Genome-wide analyses characterize shared heritability among cancers and identify novel cancer susceptibility regions.

BACKGROUND: The shared inherited genetic contribution to risk of different cancers is not fully known. In this study, we leverage results from 12 cancer genome-wide association studies (GWAS) to quantify pairwise genome-wide genetic correlations across cancers and identify novel cancer susceptibility loci. METHODS: We collected GWAS summary statistics for 12 solid cancers based on 376 759 participants with cancer and 532 864 participants without cancer of European ancestry. The included cancer types were breast, colorectal, endometrial, esophageal, glioma, head and neck, lung, melanoma, ovarian, pancreatic, prostate, and renal cancers. We conducted cross-cancer GWAS and transcriptome-wide association studies to discover novel cancer susceptibility loci. Finally, we assessed the extent of variant-specific pleiotropy among cancers at known and newly identified cancer susceptibility loci. RESULTS: We observed widespread but modest genome-wide genetic correlations across cancers. In cross-cancer GWAS and transcriptome-wide association studies, we identified 15 novel cancer susceptibility loci. Additionally, we identified multiple variants at 77 distinct loci with strong evidence of being associated with at least 2 cancer types by testing for pleiotropy at known cancer susceptibility loci. CONCLUSIONS: Overall, these results suggest that some genetic risk variants are shared among cancers, though much of cancer heritability is cancer-specific and thus tissue-specific. The increase in statistical power associated with larger sample sizes in cross-disease analysis allows for the identification of novel susceptibility regions. Future studies incorporating data on multiple cancer types are likely to identify additional regions associated with the risk of multiple cancer types.

Integrated Physiology of the Exocrine and Endocrine Compartments in Pancreatic Diseases: Workshop Proceedings.

ABSTRACT: The "Integrated Physiology of the Exocrine and Endocrine Compartments in Pancreatic Diseases" Workshop was a 1.5-day scientific conference at the National Institutes of Health (Bethesda, MD) that engaged clinical and basic science investigators interested in diseases of the pancreas. This report summarizes the workshop proceedings. The goal of the workshop was to forge connections and identify gaps in knowledge that could guide future research directions. Presentations were segregated into 6 major themes, including (a) Pancreas Anatomy and Physiology; (b) Diabetes in the Setting of Exocrine Disease; (c) Metabolic Influences on the Exocrine Pancreas; (d) Genetic Drivers of Pancreatic Diseases; (e) Tools for Integrated Pancreatic Analysis; and (f) Implications of Exocrine-Endocrine Crosstalk. For each theme, there were multiple presentations followed by panel discussions on specific topics relevant to each area of research; these are summarized herein. Significantly, the discussions resulted in the identification of research gaps and opportunities for the field to address. In general, it was concluded that as a pancreas research community, we must more thoughtfully integrate our current knowledge of the normal physiology as well as the disease mechanisms that underlie endocrine and exocrine disorders so that there is a better understanding of the interplay between these compartments.

Inferred expression regulator activities suggest genes mediating cardiometabolic genetic signals.

Expression QTL (eQTL) analyses have suggested many genes mediating genome-wide association study (GWAS) signals but most GWAS signals still lack compelling explanatory genes. We have leveraged an adipose-specific gene regulatory network to infer expression regulator activities and phenotypic master regulators (MRs), which were used to detect activity QTLs (aQTLs) at cardiometabolic trait GWAS loci. Regulator activities were inferred with the VIPER algorithm that integrates enrichment of expected expression changes among a regulator's target genes with confidence in their regulator-target network interactions and target overlap between different regulators (i.e., pleiotropy). Phenotypic MRs were identified as those regulators whose activities were most important in predicting their respective phenotypes using random forest modeling. While eQTLs were typically more significant than aQTLs in cis, the opposite was true among candidate MRs in trans. Several GWAS loci colocalized with MR trans-eQTLs/aQTLs in the absence of colocalized cis-QTLs. Intriguingly, at the 1p36.1 BMI GWAS locus the EPHB2 cis-aQTL was stronger than its cis-eQTL and colocalized with the GWAS signal and 35 BMI MR trans-aQTLs, suggesting the GWAS signal may be mediated by effects on EPHB2 activity and its downstream effects on a network of BMI MRs. These MR and aQTL analyses represent systems genetic methods that may be broadly applied to supplement standard eQTL analyses for suggesting molecular effects mediating GWAS signals.

A 584 bp deletion in CTRB2 inhibits chymotrypsin B2 activity and secretion and confers risk of pancreatic cancer.

Genome-wide association studies (GWASs) have discovered 20 risk loci in the human genome where germline variants associate with risk of pancreatic ductal adenocarcinoma (PDAC) in populations of European ancestry. Here, we fine-mapped one such locus on chr16q23.1 (rs72802365, p = 2.51 × 10-17, OR = 1.36, 95% CI = 1.31-1.40) and identified colocalization (PP = 0.87) with aberrant exon 5-7 CTRB2 splicing in pancreatic tissues (pGTEx = 1.40 × 10-69, ßGTEx = 1.99; pLTG = 1.02 × 10-30, ßLTG = 1.99). Imputation of a 584 bp structural variant overlapping exon 6 of CTRB2 into the GWAS datasets resulted in a highly significant association with pancreatic cancer risk (p = 2.83 × 10-16, OR = 1.36, 95% CI = 1.31-1.42), indicating that it may underlie this signal. Exon skipping attributable to the deletion (risk) allele introduces a premature stop codon in exon 7 of CTRB2, yielding a truncated chymotrypsinogen B2 protein that lacks chymotrypsin activity, is poorly secreted, and accumulates intracellularly in the endoplasmic reticulum (ER). We propose that intracellular accumulation of a nonfunctional chymotrypsinogen B2 protein leads to ER stress and pancreatic inflammation, which may explain the increased pancreatic cancer risk in carriers of CTRB2 exon 6 deletion alleles.

A pooled genome-wide association study identifies pancreatic cancer susceptibility loci on chromosome 19p12 and 19p13.3 in the full-Jewish population.

Jews are estimated to be at increased risk of pancreatic cancer compared to non-Jews, but their observed 50-80% excess risk is not explained by known non-genetic or genetic risk factors. We conducted a GWAS in a case-control sample of American Jews, largely Ashkenazi, including 406 pancreatic cancer patients and 2332 controls, identified in the dbGaP, PanScan I/II, PanC4 and GERA data sets. We then examined resulting SNPs with P < 10-7 in an expanded sample set, of 539 full- or part-Jewish pancreatic cancer patients and 4117 full- or part-Jewish controls from the same data sets. Jewish ancestries were genetically determined using seeded FastPCA. Among the full Jews, a novel genome-wide significant association was detected on chromosome 19p12 (rs66562280, per-allele OR = 1.55, 95% CI = 1.33-1.81, P = 10-7.6). A suggestive relatively independent association was detected on chromosome 19p13.3 (rs2656937, OR = 1.53, 95% CI = 1.31-1.78, P = 10-7.0). Similar associations were seen for these SNPs among the full and part Jews combined. This is the first GWAS conducted for pancreatic cancer in the increased-risk Jewish population. The SNPs rs66562280 and rs2656937 are located in introns of ZNF100-like and ARRDC5, respectively, and are known to alter regulatory motifs of genes that play integral roles in pancreatic carcinogenesis.

Genetic and Circulating Biomarker Data Improve Risk Prediction for Pancreatic Cancer in the General Population.

BACKGROUND: Pancreatic cancer is the third leading cause of cancer death in the United States, and 80% of patients present with advanced, incurable disease. Risk markers for pancreatic cancer have been characterized, but combined models are not used clinically to identify individuals at high risk for the disease. METHODS: Within a nested case-control study of 500 pancreatic cancer cases diagnosed after blood collection and 1,091 matched controls enrolled in four U.S. prospective cohorts, we characterized absolute risk models that included clinical factors (e.g., body mass index, history of diabetes), germline genetic polymorphisms, and circulating biomarkers. RESULTS: Model discrimination showed an area under ROC curve of 0.62 via cross-validation. Our final integrated model identified 3.7% of men and 2.6% of women who had at least 3 times greater than average risk in the ensuing 10 years. Individuals within the top risk percentile had a 4% risk of developing pancreatic cancer by age 80 years and 2% 10-year risk at age 70 years. CONCLUSIONS: Risk models that include established clinical, genetic, and circulating factors improved disease discrimination over models using clinical factors alone. IMPACT: Absolute risk models for pancreatic cancer may help identify individuals in the general population appropriate for disease interception.

Prediagnostic Leukocyte Telomere Length and Pancreatic Cancer Survival.

BACKGROUND: Leukocyte telomere length has been associated with risk of subsequent pancreatic cancer. Few prospective studies have evaluated the association of prediagnostic leukocyte telomere length with pancreatic cancer survival. METHODS: We prospectively examined the association of prediagnostic leukocyte telomere length with overall survival (OS) time among 423 participants diagnosed with pancreatic adenocarcinoma between 1984 and 2008 within the Health Professionals Follow-up Study, Nurses' Health Study, Physicians' Health Study, and Women's Health Initiative. We measured prediagnostic leukocyte telomere length in banked blood samples using quantitative PCR. Cox proportional hazards models were used to estimate HRs for OS with adjustment for potential confounders. We also evaluated 10 SNPs at the telomerase reverse transcriptase locus. RESULTS: Shorter prediagnostic leukocyte telomere length was associated with reduced OS among patients with pancreatic cancer (P trend = 0.04). The multivariable-adjusted HR for OS comparing the lowest with highest quintiles of leukocyte telomere length was 1.39 (95% confidence interval, 1.01-1.93), corresponding to a 3-month difference in median OS time. In an analysis excluding cases with blood collected within 2 years of cancer diagnosis, the association was moderately stronger (HR, 1.55; 95% confidence interval, 1.09-2.21; comparing the lowest with highest quintiles; P trend = 0.01). No prognostic association or effect modification for the prognostic association of prediagnostic leukocyte telomere length was noted in relation to the studied SNPs. CONCLUSIONS: Prediagnostic leukocyte telomere length was associated with pancreatic cancer survival. IMPACT: Prediagnostic leukocyte telomere length can be a prognostic biomarker in pancreatic cancer.

Endogenous Gastrin Collaborates With Mutant KRAS in Pancreatic Carcinogenesis.

OBJECTIVE: The KRAS gene is the most frequently mutated gene in pancreatic cancer, and no successful anti-Ras therapy has been developed. Gastrin has been shown to stimulate pancreatic cancer in an autocrine fashion. We hypothesized that reactivation of the peptide gastrin collaborates with KRAS during pancreatic carcinogenesis. METHODS: LSL-Kras; P48-Cre (KC) mutant KRAS transgenic mice were crossed with gastrin-KO (GKO) mice to develop GKO/KC mice. Pancreata were examined for 8 months for stage of pancreatic intraepithelial neoplasia lesions, inflammation, fibrosis, gastrin peptide, and microRNA expression. Pancreatic intraepithelial neoplasias from mice were collected by laser capture microdissection and subjected to reverse-phase protein microarray, for gastrin and protein kinases associated with signal transduction. Gastrin mRNA was measured by RNAseq in human pancreatic cancer tissues and compared to that in normal pancreas. RESULTS: In the absence of gastrin, PanIN progression, inflammation, and fibrosis were significantly decreased and signal transduction was reversed to the canonical pathway with decreased KRAS. Gastrin re-expression in the PanINs was mediated by miR-27a. Gastrin mRNA expression was significantly increased in human pancreatic cancer samples compared to normal human pancreas controls. CONCLUSIONS: This study supports the mitogenic role of gastrin in activation of KRAS during pancreatic carcinogenesis.

Cell-type-specific eQTL of primary melanocytes facilitates identification of melanoma susceptibility genes.

Most expression quantitative trait locus (eQTL) studies to date have been performed in heterogeneous tissues as opposed to specific cell types. To better understand the cell-type-specific regulatory landscape of human melanocytes, which give rise to melanoma but account for <5% of typical human skin biopsies, we performed an eQTL analysis in primary melanocyte cultures from 106 newborn males. We identified 597,335 cis-eQTL SNPs prior to linkage disequilibrium (LD) pruning and 4997 eGenes (FDR < 0.05). Melanocyte eQTLs differed considerably from those identified in the 44 GTEx tissue types, including skin. Over a third of melanocyte eGenes, including key genes in melanin synthesis pathways, were unique to melanocytes compared to those of GTEx skin tissues or TCGA melanomas. The melanocyte data set also identified trans-eQTLs, including those connecting a pigmentation-associated functional SNP with four genes, likely through cis-regulation of IRF4 Melanocyte eQTLs are enriched in cis-regulatory signatures found in melanocytes as well as in melanoma-associated variants identified through genome-wide association studies. Melanocyte eQTLs also colocalized with melanoma GWAS variants in five known loci. Finally, a transcriptome-wide association study using melanocyte eQTLs uncovered four novel susceptibility loci, where imputed expression levels of five genes (ZFP90, HEBP1, MSC, CBWD1, and RP11-383H13.1) were associated with melanoma at genome-wide significant P-values. Our data highlight the utility of lineage-specific eQTL resources for annotating GWAS findings, and present a robust database for genomic research of melanoma risk and melanocyte biology.

Pancreatic cancer risk is modulated by inflammatory potential of diet and ABO genotype: a consortia-based evaluation and replication study.

Diets with high inflammatory potential are suspected to increase risk for pancreatic cancer (PC). Using pooled analyses, we examined whether this association applies to populations from different geographic regions and population subgroups with varying risks for PC, including variation in ABO blood type. Data from six case-control studies (cases, n = 2414; controls, n = 4528) in the Pancreatic Cancer Case-Control Consortium (PanC4) were analyzed, followed by replication in five nested case-control studies (cases, n = 1268; controls, n = 4215) from the Pancreatic Cancer Cohort Consortium (PanScan). Two polymorphisms in the ABO locus (rs505922 and rs8176746) were used to infer participants' blood types. Dietary questionnaire-derived nutrient/food intake was used to compute energy-adjusted dietary inflammatory index (E-DII®) scores to assess inflammatory potential of diet. Pooled odds ratios (ORs) and 95% confidence intervals (CIs) were calculated using multivariable-adjusted logistic regression. Higher E-DII scores, reflecting greater inflammatory potential of diet, were associated with increased PC risk in PanC4 [ORQ5 versus Q1=2.20, 95% confidence interval (CI) = 1.85-2.61, Ptrend < 0.0001; ORcontinuous = 1.20, 95% CI = 1.17-1.24], and PanScan (ORQ5 versus Q1 = 1.23, 95% CI = 0.92-1.66, Ptrend = 0.008; ORcontinuous = 1.09, 95% CI = 1.02-1.15). As expected, genotype-derived non-O blood type was associated with increased PC risk in both the PanC4 and PanScan studies. Stratified analyses of associations between E-DII quintiles and PC by genotype-derived ABO blood type did not show interaction by blood type (Pinteraction = 0.10 in PanC4 and Pinteraction=0.13 in PanScan). The results show that consuming a pro-inflammatory diet and carrying non-O blood type are each individually, but not interactively, associated with increased PC risk.

Publisher correction: Functional characterization of a multi-cancer risk locus on chr5p15.33 reveals regulation of TERT by ZNF148.

This corrects the article DOI: 10.1038/ncomms15034.

Transcriptional regulation by NR5A2 links differentiation and inflammation in the pancreas.

Chronic inflammation increases the risk of developing one of several types of cancer. Inflammatory responses are currently thought to be controlled by mechanisms that rely on transcriptional networks that are distinct from those involved in cell differentiation. The orphan nuclear receptor NR5A2 participates in a wide variety of processes, including cholesterol and glucose metabolism in the liver, resolution of endoplasmic reticulum stress, intestinal glucocorticoid production, pancreatic development and acinar differentiation. In genome-wide association studies, single nucleotide polymorphisms in the vicinity of NR5A2 have previously been associated with the risk of pancreatic adenocarcinoma. In mice, Nr5a2 heterozygosity sensitizes the pancreas to damage, impairs regeneration and cooperates with mutant Kras in tumour progression. Here, using a global transcriptomic analysis, we describe an epithelial-cell-autonomous basal pre-inflammatory state in the pancreas of Nr5a2+/- mice that is reminiscent of the early stages of pancreatitis-induced inflammation and is conserved in histologically normal human pancreases with reduced expression of NR5A2 mRNA. In Nr5a2+/-mice, NR5A2 undergoes a marked transcriptional switch, relocating from differentiation-specific to inflammatory genes and thereby promoting gene transcription that is dependent on the AP-1 transcription factor. Pancreatic deletion of Jun rescues the pre-inflammatory phenotype, as well as binding of NR5A2 to inflammatory gene promoters and the defective regenerative response to damage. These findings support the notion that, in the pancreas, the transcriptional networks involved in differentiation-specific functions also suppress inflammatory programmes. Under conditions of genetic or environmental constraint, these networks can be subverted to foster inflammation.

Impact of Sixteen Established Pancreatic Cancer Susceptibility Loci in American Jews.

Background: The higher risk of pancreatic cancer in Ashkenazi Jews compared with non-Jews is only partially explained by the increased frequency of BRCA1 and BRCA2 mutations in Ashkenazi Jews.Methods: We evaluated the impact of 16 established pancreatic cancer susceptibility loci in a case-control sample of American Jews, largely Ashkenazi, including 406 full-Jewish pancreatic cancer patients and 2,332 full-Jewish controls, genotyped as part of the Pancreatic Cancer Cohort and Case-Control Consortium I/II (PanScan I/II), Pancreatic Cancer Case-Control Consortium (PanC4), and Resource for Genetic Epidemiology Research on Adult Health and Aging (GERA) datasets. We compared risk in full-Jewish subjects with risk in part-Jewish; non-Jewish Southern European; and in the combined non-Jewish Eastern, Northern, Southern, and Western European (non-Jewish white European) subjects from the same datasets. Jewish ancestries were genetically identified using seeded Fast principal component analysis. Data were analyzed by unconditional logistic regression, and adjusted for age, sex, and principal components.Results: One SNP on chromosome 13q22.1 (rs9543325; OR, 1.36; 95% confidence interval, 1.16-1.58; P = 10-4.1) was significant in full-Jews. Individual ORs and minor allele frequencies were similar between Jewish and non-Jewish white European subjects. The average ORs across the 16 pancreatic cancer susceptibility loci for full-Jewish, full- plus part-Jewish, non-Jewish Southern European, and non-Jewish white European subjects were 1.25, 1.30, 1.31, and 1.26, respectively.Conclusions: The 16 pancreatic cancer susceptibility loci similarly impact Jewish and non-Jewish white European subjects, both individually and as summary odds.Impact: These 16 pancreatic cancer susceptibility loci likely do not explain the higher risk seen in Ashkenazi Jews. Cancer Epidemiol Biomarkers Prev; 26(10); 1540-8. ©2017 AACR.